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Abstract — Finite alpliabet iterative decoders (FAIDs) for LDPC 
codes were recently sliown to be capable of surpassing the 
Belief Propagation (BP) decoder in the error floor region on the 
Binary Symmetric channel (BSC). More recently, the technique 
of decimation which involves fixing the values of certain bits 
during decoding, was proposed for FAIDs in order to make 
them more amenable to analysis while maintaining their good 
performance. In this paper, we show how decimation can be used 
adaptively to further enhance the guaranteed error correction 
capability of FAIDs that are already good on a given code. 
The new adaptive decimation scheme proposed has marginally 
added complexity but can significantly improve the slope of the 
error floor performance of a particular FAID. We describe the 
adaptive decimation scheme particularly for 7-level FAIDs which 
propagate only 3-bit messages and provide numerical results for 
column-weight three codes. Analysis suggests that the failures of 
the new decoders are linked to stopping sets of the code. 

I. Introduction 

The error floor problem of low-density parity-check (LDPC) 
codes under iterative decoding is now a well-known problem, 
where the codes suffer from an abrupt degradation in their 
error-rate performance in spite of having good minimum 
distance. The problem has been attributed to the presence 
of harmful configurations generically termed as trapping sets 
|[T] present in the Tanner graph, which cause the iterative 
decoder to fail for some low-noise configurations, thereby 
reducing its guaranteed error correction capability to an extent 
that is far from the limits of maximum likelihood decoding. 
More importantly for the BSC, the slope of the error floor is 
governed by the guaranteed correction capability |2|. 

Recently, a new class of finite alphabet iterative decoders 
(FAIDs) that have a much lower complexity than the BP 
decoder, were proposed for LDPC codes on the BSC lO lH 
and were shown to be capable of outperforming BP in the 
error floor Numerical results on several column-weight three 
codes showed that there exist 7-level FAIDs requiring only 
3 bits of precision that can achieve a better guaranteed error 
correction ability than BP, thereby surpassing it in the error 
floor region. However, analyzing these decoders for providing 
performance guarantees still proved to be difficult. 

More recently, decimation-enhanced FAIDs |5| were pro- 
posed for BSC in order to make FAIDs more amenable to anal- 
ysis while maintaining their good performance. The technique 
of decimation involves guessing the values of certain bits, and 



fixing them to these values while continuing to estimate the 
remaining bits (see Q for references). In fSl, the decimation 
was carried out by the FAID based on messages passed for 
some iterations, and a decimation scheme was provided such 
that a 7-level DFAID matched the good performance of the 
original 7-level FAID while being analyzable at the same time. 

In this paper, we show how decimation can be used 
adaptively to further increase the guaranteed error correction 
capability of FAIDs. The adaptive scheme has only marginally 
increased complexity, but can significantly improve the error- 
rate performance compared to the FAIDs. We specifically 
focus on decoders that propagate only 3-bit messages and 
column-weight three codes since these enable simple imple- 
mentations and thus have high practical value. We also provide 
some analysis of the decoders which suggests that the failures 
are linked to stopping sets of the code. Numerical results are 
also provided to validate the efficacy of the proposed scheme. 

II. Preliminaries 

Let G = {V U C,E) denote the Tanner graph of an {n,m) 
binary LDPC code C with the set of variable nodes V = 
{vi, ■ ■ ■ ,Vn} and set of check nodes C — {ci, • • • ,Cm}- E is 
the set of edges in G. A code C is said to be dy-left-regular if 
all variable nodes in V of graph G have the same degree dy. 
The degree of a node is the number of its neighbors. c?„ii„ is 
the minimum distance of the code C. 

A trapping set is a non-empty set of variable nodes in G 
that are not eventually corrected by the decoder HI. 

A multilevel FAID ^ is a 4-tuple given by ^ = 
{Ai,y,^v,^c) \M- The messages are levels confined to an 
alphabet M of size 2s + 1 defined as = {0, ±Li : I < i < 
s}, where Li £ M+ and Li > Lj for any i > j. y denotes the 
set of possible channel values that are input to the decoder 
For the case of BSC as 3^ = {±C : C G M+}, and for each 
variable node Vi E V, the channel value yt E y is determined 
by Vi — (~1)'^*C, where is received from the BSC at Vi. 

Let TOi, • • • , m„ denote the incoming messages to a node 
with degree u. $c : Ai'^"^^ is the update function used 

at a check node with degree dc and defined as 

^c{mi,...,md^_i) = \ TT sgn(TOj) min i\mj\) 

\ , / i6{l,...,dc-l} 
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where sgn denotes the standard signum function. 

$1, : y X Ai'^^'^^ — > is a symmetric update rule used at 
a variable node with degree dy and is defined as 

^viyi,rni,m2, ■ ■ ■ ,TOd„-i) = Q ^ rrij + ojc ■ Vt \ ■ 

The function Q is defined based on a threshold set 3^ ~ {Ti : 
1 < i < 3 + 1} where T, e M+ and T, > Tj if i > j, and 
Ts+i — oo, such that Q{x) = sgn{x)Li if Ti < \x\ < T^+i, 
and Q{x) = otherwise. The weight uic is computed using 
a symmetric function ft : J\4'^^~'^ — >■ {0, 1}. Based on this, 
$^ can be described as a linear-threshold (LT) or non-linear 
threshold (NLT) function. If f2 = 1 (or constant), then it is 
an LT function, else it is an NLT function. $„ can also be 
described as a look-up table (LUT) as shown in Table |l] (for 
Ui = — C, it can be obtained from symmetry). This rule will 
be one of the rules used in the proposed decimation scheme. 

Let T/'IG) denote the computation tree of graph G cor- 
responding to a decoder ^ enumerated for k iterations with 
variable node Vi £ V as its root. A node w E T^{G) is a 
descendant of a node u E T^{G) if there exists a path starting 
from node w to the root Vi that traverses through node u. 

Definition 1 (Isolation assumption): Let _ff be a subgraph 
of G induced hy P QV with check node set W C C. The 
computation tree T^{G) with the root Uj G P is said to be 
isolated if and only if for any node m ^ P U in T/'iG), u 
does not have any descendant belonging to P U W . If TJ*^ (G) 
is isolated Vu^ G P, then the subgraph H is said to satisfy the 
isolation assumption in G for k iterations. 

Remark: The above definition is a revised version of the 
one given in f3l. 

The critical number of a FAID on a subgraph H is the 
smallest number of errors for which F fails on H under the 
isolation assumption. 

Let M{u) denote the set of neighbors of a node u in the 
graph G and let M{U) denote the set of neighbors of all u e 
U. Let mk{vi,Af{vi)) denote the set of outgoing messages 
from Vi to all its neighbors in the A;*'' iteration. Let bi denote 
the bit associated to a variable node Vi E V that is decided by 
the iterative decoder at the end of the k*^ iteration. 

III. Adaptive decimation-enhanced faids 

We will first provide some definitions and notations related 
to the concept of decimation and discuss its incorporation into 



the framework of FAIDs before we delve into the details of 
adaptive decimation. 

Definition 2: A variable node Vi is said to be decimated 
at the end of iteration if 6f is set to b* Vfc > I. Also, 
mk{vi,Af{vi)) = Lg}, Vfc > / irrespective of its 

incoming messages, i.e., Vi always sends the strongest possible 
messages. 

Remark: If a node Vi is decimated, then all its descendants 
in the computation tree Tl'{G) can be deleted since the node 
always sends (—!)''» to its parent. 

A decimation rule /3 : y x A^'^" — > {—1, 0, 1} is a function 
used by the decoder to decide whether a variable node should 
be decimated and what value it should be decimated to. Let 
7i denote the output of a decimation rule applied to a node 
Vi. If 7i = 0, then the node is not decimated. If 7^ — 1, then 
b* = 0, and if 7i = -1, then b* = 1 

Remark: The decimation rule is a function of the channel 
value and the most recent incoming messages received by the 
variable node before the decimation rule is applied. 

We shall refer to each instance of applying a decimation 
rule on all the variable nodes as a decimation round. 

There are two key aspects to note regarding the application 
of a decimation rule: 1) the decimation rule is applied after 
messages are passed iteratively for some I iterations, and 
2) after each instance of applying the decimation rule, all 
messages are cleared to zero (which is practically restarting the 
decoder except that the decimated nodes remain decimated). 

Let denote the number of decimation rounds carried out 
by the decoder with a given decimation rule /3 beyond which 
no more variable nodes are decimated. 

Definition 3: The residual graph G' is the induced sub- 
graph of the set of variable nodes in G that are not decimated 
after Nd decimation rounds. 

We can now formally define the class of adaptive 
decimation-enhanced multilevel FAIDs (ADFAIDs) as follows. 
A decoder belonging to such a class denoted by is defined 
as = (7W,3^,<I>^,$$5,$;,6,<I>c), where the sets M and 
y, and the map $c are same as the ones defined for a multi- 
level FAID. The map $^ : 3^ x M'^^'^ x {-1, 0,1} M is 
the update rule used at the variable node. It requires the output 
of a decimation rule /3 as an one of its arguments and also uses 
the maps $^ : 3^ x ^ X and : 3^ x M'^^'^ M 

to compute its output. For simplicity, we define it for the case 
of = 3 as follows. 

r $^(2/„mi,TO2), if 7i=0,P< Nd 

I liLs, if 7i = ±1 

where p denotes the p*'* decimation round completed by the 
decoder The maps <i>^ and $Jj are defined as either LT or NLT 
functions or as look-up tables similar to of a FAID 

Remark: The new class of decoders proposed in this paper 
use two different maps, ^f, and for updating the messages 
on non-decimated variable nodes. is the map used to 
update messages specifically during the decimation procedure. 



whereas $[, is the map used to decode the remaining non- 
decimated nodes after the decimation procedure is completed. 
Also note that for the case of |A^| = 7, we restrict the defini- 
tion of <S>^ to satisfy $^(C,0,0) = Li, <^>^{C, Li, Li) = La, 
and <I>55(C, L2, -^2) = ^3- is also defined similarly. 

Proposition 1: Given a decimation rule (3, if the number of 
decimated nodes after the p*'* decimation round is the same as 
the number of decimated nodes after the {p— 1)*'' decimation 
round, then no additional nodes will get decimated in the 
subsequent decimation rounds. 

Remark: This would be the stopping criterion used for the 
decimation procedure. In the above case, iV^ — p. 

The set B is the set of decimation rules used for adaptive 
decimation and for any /3 e S, it satisfies the following 
properties (specified for — 3). 

1) /3(C,mi,m2,m3) = -^(-C, -toi, -ms, -mg) 
Vmi, m2, ma € M 

2) ^(C,TOi,m2,m3) 7^ -1 and /3(-C, mi, 777,2, ma) 7^ 1 
Vmi, m2, ma e 

3) Given mi, m2, ma G 7M, if /3(C, mi, TO2, ma) = 1, then 
/3(C, m'l, mj, mg) = 1 Vm'j,m2,m3 S such that 
m'l > mi, 771-2 — ™2, and mg > mg. 

Remark: Property 2 implies that a node 7;^ can be decimated 
to zero only if 7/,; = C and to one only if 7/i = — C. 
Consequently a node initially correct will never be decimated 
to a wrong value and a node initially wrong will never be 
decimated to the correct value. Then, the necessary condition 
for successful decoding is that no node initially in error is 
decimated. We shall now restrict our discussion to — 3 for 
the remainder of the paper. 

For a given decimation rule /3, a set S can be used to com- 
pletely specify /3, where S is defined as the set of all unordered 
triples (mi,m2,ma) G M.^ such that /3(C, mi, m2, ma) = 
1. Note that for any unordered triple (mi,m2,ma) G S, 
/3(— C, — mi, — m2, — ma) = —1 by property 1, so S is 
sufficient to completely specify /3. A /? is considered to be a 
conservative decimation rule if |S| is small and an aggressive 
rule if |S| is large. 

Note that the class of decimation-enhanced FAIDs defined 
in our previous work ||5l is a special case of the newly 
proposed decoders where $j5 = <i>J^ and B — {/?}. In other 
words, only a single non-adaptive decimation rule and a single 
map is used for updating messages in the DFAIDs of Q. 

For the remainder of the paper, we shall refer to variable 
nodes that are initially in error in G as error nodes and variable 
nodes that are initially correct as correct nodes. 

A. Motivation for adaptive decimation 

Given an error pattern of relatively low weight (< 
l^ rfmin-i j-)^ the primary role of decimation is to isolate the 
subgraph associated with the error pattern from the rest 
of the graph by decimating as many correct nodes outside 
this subgraph as possible. The rationale behind resetting the 
messages to zero at the end of each decimation round is 
to allow more non-decimated correct nodes that are close to 
the neighborhood of the decimated correct nodes to possibly 



be decimated as long as none of the error nodes have been 
decimated. This is possible since the decimated nodes always 
send the strongest message (±Ls)- 

Now if a given error pattern is such that the error nodes 
are relatively clustered with many interconnections between 
them through their neighboring check nodes, then a more 
conservative /3 would have to be used by the decoder to ensure 
that none of the error nodes are decimated. However, if the 
error pattern is such that the error nodes are more spread 
out, then it may be desirable to use a more aggressive j3 as 
there will be many correct nodes in the neighborhood of the 
error nodes that can be decimated without decimating the error 
nodes, and, in turn, possibly help the decoder to converge. This 
is our main motivation for the use of adaptive decimation in 
the newly proposed decoders, and we will eventually show 
that adaptive decimation can help achieve an increase in the 
guaranteed error correction capability of the code. 

B. Proposed scheme 

We will now describe a particular adaptive decimation 
scheme used by the decoder in order to enhance the 
guaranteed error correction capability. In the proposed scheme, 
the set B consists of two decimation rules, namely B = 
{/?^^\ /3^^'}, where S^^' and 'EP'^ are the sets of unordered 
triples that completely specify the rules /J^^' and (3^"^^ respec- 
tively. The rule /J^^' is used only once at the end of the third 
iteration, and then after that, /J^^^ is used after every two 
iterations {I = 2). The use of adaptive decimation is carried 
out only through as follows. 

We define a sequence of decimation rules from 
by considering subsets of S*^^-' with increasing size. 
Let be the number of rules in the sequence {Z?^^-''-''},- 
and let S^^'I^l denote the set that specifies the rule /S^^^I^l. 
Then S^^^I^l is defined for each /^t^)!-?! in a way such that 

s(2)b] c V7 e {1, . . . , with = 

This implies that the sequence of rules are such that 
is less conservative than with being the most 

conservative and /3(2)[A'/3] — ^(2) being least conservative (or 
most aggressive). Note that the each subset S'^^^I-'l must be 
chosen in a manner that ensures that its corresponding rule 
^(2)b] satisfies the properties of j3 mentioned previously. 

For a given error pattern, the decoder starts the decimation 
procedure by passing messages using the map $^ and applying 
the decimation rule /3'^' at the end of the third iteration 
after which the messages are reset to zero. Then the most 
conservative rule in the sequence which is 

is used after every two iterations (followed by resetting the 
messages) until no more nodes can be decimated. The map $JJ 
then is used to decode the remaining non-decimated nodes. If 
the decoder still does not converge, then the whole decoding 
process is repeated by using a more aggressive rule /3(^)[^1 
in place of Z?*^^^!^!. This decoding process continues until the 
decoder converges or until all rules in the sequence {/3^^'''-''}j 
have been used. Let Ni, denote the number of decimated bits 
at the end of a decimation round. The decoding scheme can 



be summarized as follows. Note that this scheme is devised 
particularly for the case of = 7. 

Algorithm 1 Adaptive decimation-enhanced FAID algorithm 

1) Set j = 1. Note that <i>c will always be used to update 
messages at the check node. 

2) Initialize ji — Vvi S V. 

3) Start the decimation procedure by passing messages for 
three iterations using If the decoder converges within 
those three iterations, STOP. 

4) Apply decimation rule (3^-^^ for every Vi e V. Then reset 
all messages to zero and set q — 0. 

5) Pass messages for two iterations using ^f, for update at 
the non-decimated nodes. If the decoder converges within 
those two iterations, STOP. 

6) Apply decimation rule /S^^^I^l only on nodes Vi for which 
7i = 0. Then reset all messages to zero. If Ni, > q, 
q ~ Nb and go back to step 5, else go to step 7. 

7) Pass messages using $^ on the nodes Vi for which 7; = 0. 

8) If decoder converges or has reached maximum allowed 
iterations, STOP. Else j = j + 1. 

9) If j > Np STOP Else go to step 2. 



Remarks: 1) The only criteria used by the decoder to decide 
when to use a more aggressive rule 

on a given error 

pattern is whether the decoding has failed. 2) The reason for 
applying p'^^^ at the end of third iteration is that at least three 
iterations are required for a ±£3 to be passed. 3)The reason for 
the choice of every 2 iterations for applying /^'^^^[■'l is because 
it is small which helps prevent the growth of wrong message 
strengths but just enough to allow all levels in M. to be passed. 

C. Choice of $^ and $^ 

For the proposed decoders, the map $Jj is simply chosen to 
be the of a particular FAID already known to be good on a 
given code, and for which we want to improve the guaranteed 
error correction capability. For the numerical results, <I>Jj is 
chosen to be the <I>^ of a 7-level FAID defined by Table |l] 

The choice of on the other hand is non-trivial. It is 
designed based on analyzing messages that are passed within 
dense subgraphs that could potentially be trapping sets for 
a given FAID when errors are introduced in them under the 
isolation assumption. The rule is chosen under the premise that 
the growth of message strengths within the subgraph should 
be slow since many correct nodes in the subgraph would most 
likely be connected to error nodes, and multiple error nodes 
may be interconnected to each other in the subgraph (if the 
number of errors introduced is comparable to the size of the 
subgraph). Explicit design methods for $^ are not discussed in 
this paper, but we provide a particular that was designed 
based on the above philosophy and used for the numerical 
results. It is an LT rule (see Section [ll|, so it can be described 
by assigning values to thresholds in and messages in M.. 
The map is defined with the following assignments; Li = 1.1, 
L2 = 2.3, L3 = 6.6, Ti = 0.8, T2 = 2.8, T3 ^ 4. This was 
found to be a good rule for decimation. 



D. Analysis 

Due to page constraints, no proofs are provided (but they 
will be provided in the journal version of this paper). For the 
analysis, we assume that the all-zero codeword is transmitted 
which is valid since the decoders considered are symmetric. 

Proposition 2: A node Vi can receive a ±£3 from its 
neighbor Cj in the first or second iteration after resetting the 
messages, only if all nodes in Af{cj)\vi have been decimated. 

Lemma 1: If /3(2)M (C, L3, -ij, -L2) = 1 and if Vc^ G 
N{vi) all error nodes in Af{ck)\vi are non-decimated, then a 
correct node Vi will be decimated if it receives an L3 during 
a decimation round. 

Remark: Note that /3*^'^'[-'l will always be defined so that 
^(2)[i] _272, — L2) = 1 for any j as explained in the 
next subsection. Also note how resetting messages at the end 
of each decimation round can help with decimating more 
correct nodes due to the above lemma. 

Theorem 1: If (C, £3, -L2, -^2) = 1 and no error 

node is decimated, then any correct node in the residual graph 
G' is connected to check nodes that have at least degree-two. 

Corollary 1: If Theorem 1 holds and no error node in the 
residual graph G' is connected to a degree-one check node, 
then G' is a stopping set. 

Remark: Note that if an error node in the residual graph 
G' is connected to a degree-one check node, it would receive 
Lj, in every iteration for the remainder of the decoding (again 
assuming no error nodes are decimated), and this will most 
likely lead to a decoder convergence. Therefore, if no error 
node is decimated, the decoder is more likely to fail when the 
residual graph G' is a stopping set (refer to fT'l for def.). 

The above remark is an important observation since we can 
now design the rules /S^^' and the sequence {/3^^'""''}j based 
on analyzing error patterns whose errors are entirely contained 
in the minimal stopping sets of a given code. For instance, if 
our goal is to correct up to t-errors, then we consider all error 
patterns up to a weight t in the stopping sets in order to design 
and 

If FAID with <!>„ = $^ has a critical number of t+1 
on a stopping set whose induced subgraph is iJ, then is 
guaranteed to correct up to t errors introduced in H on the 
code if the residual graph is H. In other words, if "I>'^ has a 
critical number of t+1 on a stopping set, it is more likely to 
correct up to t errors on the code containing the stopping set. 

E. Discussion on design of decimation rules and /J^^-* 

The design of /J^^^ involves selecting the triples that should 
be included in 2*^^^ which depends on the number of errors 
we are trying to correct and the type of harmful subgraphs 
present in G. should be chosen to be conservative enough 
so that no error nodes are decimated. On the other hand, 
the design of /^^^^^ not only involves selecting the triples that 
should be included in 'E!^'^\ but also determining a specific 
ordering on the triples that will be included in subsets S'^'^'I^l 
which determine the sequence of rules used starting 

from the least conservative rule, and this is dependent on the 
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Fig. 1. FER performance comparison on the (155,93) Tanner code 



Structure of the code. Both rules can be designed by analyzing 
them on errors introduced in stopping sets of the code. 

In order to specify the set S^^^^ we just specify the message 
triples with the weakest values. For specifying S^^^ in a 
concise way, we shall introduce some notations. Let S*^^^ be 
divided into two disjoint subsets, i.e., S'^'^' = A U F, where A 
is a subset that contains all triples {L^, 1712,1713) e such 
that 1712,1713 > —L2- Based on the analysis in the previous 
subsection, any S^^^I^l defined should always have A as its 
subset, regardless of the code. 

IV. Numerical Results and Discussion 

Numerical results are provided in Fig. [T] and Fig. |2] for two 
codes: the well-known (155,93) Tanner code and a structured 
rate 0.75 (732,183) code constructed based on latin squares ||6l 
with dmin = 12. For the Tanner code, the set S*^^^ contains 
all triples (mi, TO2, ms) G Ai^ such that (7711,7712,^3) > 
(L3,0,0) and (7771,7712,7773) > {L2,L2,Li) (comparison 
is componentwise). For the high-rate structured code, S'^-' 
contains all triples such that (7771,7712,^3) > {L^, Li, —L^), 
(7771,7772,7773) > (L3,-Li,-Li), and (7771,7772,7773) > 
{L2,Li,Li). = 12 for the Tanner code and = 24 

for the (732,183) code. For the set S^^) = A U T, the 
set F that was used for the structured code and the Tanner 
code is given in Tables [ll] and III resp. The cardinalities 
of the subsets of S^^^ used by each of the two codes are 
= {24,25,26,27,28,29,30,31,32,33,34,35} 
and {|S(2)W|}^- = {23,25,26,27,29} resp. The maximum 
number of iterations allowed for BP and 7-level FAID, and 
for of the 7-level ADFAID, was 100. 

The significant improvement in the slope of the error floor 
by using the 7-level ADFAID is evident. For the Tanner code, 
it was verified that all 6-error patterns are corrected by the 
7-level ADFAID while the 7-level FAID corrects all 5 -errors 
and BP fails on 5-errors. For the high-rate structured code, 
no failed 5 -error patterns were found either in the region of 
simulation shown in Fig. [2] or by a brute force check (still 
running), which is significant since the code has dmin — 12. 
This shows that for certain high-rate codes whose graphs are 
relatively dense and for which it becomes difficult to ensure 
high dmin in the code, the FAIDs with adaptive decimation 
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can possibly come close to achieving the guaranteed error 
correction of maximum likelihood decoding. Note that the 7- 
level ADFAIDs are still 3-bit message passing decoders which 
have reasonable complexity, and that is still lower than BP. 

Acknowledgment 

This work was funded by the NSF under the grants CCF- 
0830245 and CCF-0963726, and by the grant "Institut Uni- 
versitaire de France". 



References 



[1] 



T. Richardson, "Error floors of LDPC codes," in Proc. 41st Annual 
AUerton Conf. on Commun., Control and Computing, 2003. 
[2] M. Ivkovic, S. K. Chilappagari, B. Vasic, "Eliminating trapping sets in 
low-density parity-check codes by using Tanner graph covers," IEEE 
Trans, on Inf. Theory, vol. 54, no. 8, pp. 3763-3768, 2008. 
S. K. Planjery, D. Declercq, S. K. Chilappagari, B. Vasic, "Multilevel de- 
coders surpassing belief propagation on the binary symmetric channel," 
Proc. Int. Symp. on Inform. Theory,pp. 769-773, Jul. 2010. 
S. K. Planjery, D. Declercq, L. Danjean, B. Vasic, "Finite alphabet 
iterative decoders for LDPC codes surpassing floating-point iterative 
decoders," Electron. Lett. vol. 47, no. 16, pp. 919-921, Aug. 2011. 
S. K. Planjery,B. Vasic,D. Declercq, "Decimation-enhanced finite alpha- 
bet iterative decoders for LDPC codes on the BSC," Proc. Int. Symp. Inf. 
Theory (ISIT20I1}, pp. 2383-2387, St. Petersburg, Russia, Jul. 2011. 
D. V. Nguyen, B. Vasic, M. Marcellin, S.K. Chilappagari,"On the 
construction of structured LDPC codes free of small trapping sets," IEEE 
Trans. Inf. Theory {accepted for publication) 

A. Orlitsky, K. Vishwanathan, J. Zhang,"Stopping set distribution of 
LDPC code ensembles," IEEE Trans. Inf. Theory, vol. 51, pp. 929-953, 
Mar. 2005. 



[3] 



[4] 



[5] 



[6] 



[7] 



